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Abstract  Sulfur  mustard  [SM,  bis-(2-chloroethyl) 
sulfide]  is  a  highly  reactive  bifunctional  alkylating  agent 
that  has  been  used  as  a  vesicating  agent  in  warfare  scenar¬ 
ios  to  induce  severe  lung,  skin,  and  eye  injury.  SM  cuta¬ 
neous  lesions  are  characterized  by  both  vesication  and 
severe  inflammation,  but  the  molecular  mechanisms  that 
lead  to  these  signs  and  symptoms  are  not  well  understood. 
There  is  a  pressing  need  for  effective  therapeutics  to  treat 
this  injury.  An  understanding  of  the  molecular  mecha¬ 
nisms  of  injury  and  identification  of  potential  therapeutic 
targets  is  necessary  for  rational  therapeutic  development. 
We  have  applied  a  high-throughput  small  interfering  RNA 
(siRNA)  screening  approach  to  the  problem  of  SM  cuta¬ 
neous  injury  in  an  effort  to  meet  these  needs.  Our  siRNA 
screening  efforts  have  initially  focused  on  SM-induced 
inflammation  in  cutaneous  injury  since  chronic  inflamma¬ 
tion  after  exposure  appears  to  play  a  role  in  progressive 
clinical  pathology,  and  intervention  may  improve  clinical 
outcome.  Also,  targets  that  mitigate  cellular  injury  should 
reduce  the  inflammatory  response.  Historical  microarray 
data  on  this  injury  were  mined  for  targets  and  pathways 
implicated  in  inflammation,  and  a  siRNA  library  of  2,017 
targets  was  assembled  for  screening.  Primary  screening 
and  library  deconvolution  were  performed  using  human 
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HaCaT  keratinocytes  and  focused  on  cell  death  and 
inflammatory  markers  as  end  points.  Using  this  in  vitro 
approach,  we  have  identified  and  validated  novel  targets 
for  the  potential  treatment  of  SM-induced  cutaneous 
injury. 

Keywords  Sulfur  mustard  •  Cutaneous  injury  •  siRNA  • 
High-throughput  screening  •  Inflammation  •  In  vitro  models 

Introduction 

Sulfur  mustard  (SM,  bis-(2-chloroethyl)  sulfide)  is  a 
highly  reactive  bifunctional  alkylating  agent  and  has  been 
used  in  warfare  scenarios  as  a  potent  vesicating  agent  to 
induce  severe  lung,  skin,  and  eye  injury  (for  a  review,  see 
Papirmeister  et  al.  1991).  SM  covalently  modifies  DNA, 
proteins,  and  other  macromolecules  (Deaton  et  al.  1990; 
Dillman  et  al.  2003;  Langenberg  et  al.  1998;  Papirmeister 
et  al.  1991).  DNA  damage  has  long  been  hypothesized  to 
play  a  primary  role  in  the  cellular  response  to  SM  exposure 
(Fox  and  Scott  1980;  Herriott  and  Price  1948;  Lawley  and 
Brookes  1965,  1968;  Papirmeister  et  al.  1969);  however, 
the  contribution  of  other  alkylation  events  to  the  response 
remains  largely  unknown.  Clinically,  SM  cutaneous  lesions 
are  characterized  by  both  vesication  and  severe  inflamma¬ 
tion  (Papirmeister  et  al.  1991).  The  molecular  mechanisms 
that  lead  to  these  signs  and  symptoms  are  not  well  under¬ 
stood.  Microarray  studies  have  shown  that  several  molec¬ 
ular  pathways  are  activated  in  response  to  SM  exposure 
(Dillman  et  al.  2006;  Rogers  et  al.  2004;  Sabourin  et  al. 
2004).  These  can  be  categorized  broadly  into  inflammatory 
and  apoptotic/necrotic  pathways.  Several  other  studies  have 
implicated  a  number  of  different  pathways  and  signaling 
molecules  such  as  p53,  p38a  MAPK,  and  caspases  just  to 
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name  a  few  (for  a  review,  see  Kehe  et  al.  2009;  Ruff  and 
Dillman  2007).  However,  inhibitors  of  these  targets  have 
not  been  very  successful  in  ameliorating  the  effects  of  SM 
exposure. 

The  molecular  and  cellular  response  to  any  disease  or 
injury  is  typically  complex;  however,  in  many  cases,  these 
complexities  often  follow  a  relatively  defined  initiating 
mechanism.  In  the  case  of  SM  injury,  even  the  initiat¬ 
ing  events  are  complex  because  SM  reacts  with  numer¬ 
ous  biological  molecules.  This  myriad  of  other  molecules 
alkylated  by  SM  likely  contributes  to  the  clinical  progres¬ 
sion  of  injury  despite  the  long-held  hypothesis  that  SM 
DNA  damage  is  the  primary  lesion  leading  to  cell  death. 
The  regulatory  pathways  involved  in  response  to  SM  are 
also  complex  and  show  interplay  between  major  regula¬ 
tors  of  major  biological  pathways  such  as  inflammation, 
cellular  differentiation,  and  cell  death.  For  example,  using 
siRNA,  our  research  group  identified  the  MAPK  p38a 
signaling  pathway  as  a  critical  regulator  of  vesicant- 
induced  inflammatory  cytokine  production  (Ruff  and 
Dillman  2010).  However,  p38  has  also  been  shown  to  pro¬ 
mote  SM-induced  keratinocyte  terminal  differentiation, 
a  role  that  is  negatively  influenced  by  ERK1/2  activity 
(Popp  et  al.  2011).  There  is  a  pressing  need  to  understand 
the  molecular  mechanisms  of  this  injury  and  to  identify 
potential  therapeutic  targets.  Given  the  diversity  of  mol¬ 
ecules  that  react  with  SM  and  the  complexity  of  the  regu¬ 
latory  pathways  involved,  traditional  methods  of  pathway 
mapping  are  unlikely  to  meet  these  needs  in  a  timely 
manner. 

Our  earlier  work  on  the  role  of  p38a  in  SM-induced 
inflammation  demonstrated  that  siRNA  approaches  can 
elucidate  mechanisms  of  vesicant  toxicity  and  poten¬ 
tially  identify  targets  for  drug  screening.  High-throughput 
siRNA  screening  is  used  by  the  pharmaceutical  industry 
to  mechanistically  assess  gene  function  and  identify  ther¬ 
apeutic  targets.  For  these  reasons,  we  have  applied  high- 
throughput  siRNA  screening  technology  to  the  problem  of 
SM  cutaneous  injury.  Our  siRNA  screening  efforts  have 
initially  focused  on  SM-induced  inflammation  in  cutane¬ 
ous  injury.  Historical  microarray  data  on  this  injury  were 
mined  for  targets  and  pathways  implicated  in  inflamma¬ 
tion.  From  these  data,  a  siRNA  library  of  2,017  targets 
was  assembled  for  screening.  Primary  screening  and 
library  deconvolution  were  performed  using  HaCaT  cells, 
a  p5 3 -mutated  immortal  human  keratinocyte  line  (Bouk- 
amp  et  al.  1988).  Cells  were  transfected  with  siRNA  prior 
to  exposure  to  SM,  and  the  effects  on  cytokine  produc¬ 
tion  and  cell  viability  were  analyzed.  Using  this  in  vitro 
approach,  we  have  identified  and  validated  novel  targets 
that  regulate  SM-induced  inflammation  that  may  poten¬ 
tially  be  useful  for  the  treatment  of  SM-induced  cutane¬ 
ous  inflammation. 


Materials  and  methods 

Cell  culture 

HaCaT  cells,  a  p5 3 -mutated  immortal  human  keratinocyte 
line,  were  a  kind  gift  from  Dr.  John  Schlager.  These  cells 
were  selected  for  study  because  they  are  a  human  model 
and  have  been  historically  used  to  study  SM  in  vitro. 
HaCaT  cells  were  grown  in  Dulbecco’s  modified  Eagle 
medium  (DMEM)  (Life  Technologies,  Grand  Island,  NY) 
supplemented  with  10  %  fetal  bovine  serum  and  1  %  L-glu- 
tamine  (Life  Technologies).  Cells  were  incubated  at  37  °C 
with  5  %  C02,  maintained  in  7 5 -cm2  cell  culture  flasks, 
and  grown  to  70-80  %  confluency  before  trypsin  detach¬ 
ment  and  reseeding  at  2.5  x  103  cells/cm2.  For  primary 
screening  and  library  deconvolution,  cells  were  trypsinized 
and  seeded  into  96- well  plates  at  3,000  cells/well  using  a 
HydraDT  96  channel  pipetting  robot  fitted  with  an  XYZ 
stage  (ThermoFisher,  Chicago,  IL).  For  siRNA  validation, 
cells  were  seeded  into  6- well  plates  at  50,000  cells/well 
using  a  standard  laboratory  pipette. 

siRNA  and  transfection 

Dharmacon  siGENOME  SmartPool  siRNA  was  used  for  pri¬ 
mary  screening  and  Dharmacon  ON-TARGETplus  siRNA 
was  used  for  library  deconvolution  and  target  validations  that 
utilized  siRNA.  Lipofectamine  RNAiMAX  transfection  rea¬ 
gent  (Life  Technologies)  was  used  to  transfect  the  siRNA.  For 
screening  and  deconvolution,  the  siRNA  libraries  were  recon¬ 
stituted,  and  transfection  mixes  were  prepared  and  added  to 
cells  using  a  ThermoFisher  HydraDT  96  channel  pipetting 
robot.  For  target  validation,  siRNA  was  purchased  in  tube 
format.  After  reconstitution  of  the  siRNA,  transfection  mixes 
were  prepared  and  added  to  cells  using  standard  single  chan¬ 
nel  digital  laboratory  pipettes.  All  siRNA  was  reconstituted 
with  Thermo  Scientific  5X  siRNA  buffer.  siRNA  and  trans¬ 
fection  reagent  dilutions  for  transfection  were  prepared  with 
Opti-MEM-reduced  serum  medium  (Life  Technologies). 
Transfections  were  performed  1  day  after  cell  seeding.  Cells 
were  transfected  for  6  h,  washed  with  phosphate-buffered 
saline  (PBS,  10  mM  sodium  phosphate,  0.9  %  NaCl,  pH 
7.4),  and  overlaid  with  fresh  DMEM  culture  medium  supple¬ 
mented  with  1  %  Penicillin/Streptomycin  solution  (Sciencell, 
Carlsbad,  CA).  Cells  were  fed  the  next  day  with  antibiotic- 
free  DMEM,  and  exposures  were  performed  the  following 
day.  This  provided  a  two-day  post-transfection  time  period 
prior  to  experimentation  to  allow  for  maximum  target  knock¬ 
down.  For  primary  screening  and  deconvolution,  each  siRNA 
was  transfected  in  six  replicate  wells  on  six  different  plates. 
Three  replicate  plates  were  used  for  SM  exposures,  and  three 
were  unexposed  controls.  Controls  on  each  96- well  plate 
included  non-targeting  siRNA,  two  toxic  siRNAs  [siTOX  and 
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siRNA  targeting  Kif  11  (Dharmacon),  positive  control  p38a 
Validated  Stealth  RNAi  (Life  Technologies)],  OptiMEM 
transfection  diluent,  and  naive  cells.  Because  siRNA  valida¬ 
tion  studies  were  performed  in  6-well  plates,  controls  were  on 
plates  separate  from  the  experimental  siRNAs  and  included 
non-targeting  siRNA,  positive  control  p38a  Validated  Stealth 
RNAi,  and  naive  cells. 

In  situ  mRNA  quantitation  and  high-content  analysis 

Transfection  conditions  for  individual  targets  were  opti¬ 
mized  following  deconvolution  of  the  library  using  the 
Quantigene  View  assay  (Affymetrix,  Santa  Clara,  CA)  to 
assess  gene  expression  levels.  The  assay  was  performed 
according  to  the  manufacturer’s  instructions.  Briefly,  cells 
were  fixed  using  a  2  %  formaldehyde  solution.  Cells  were 
then  permeabilized  and  digested  with  the  included  protease 
solution  diluted  1:1,000  before  hybridization  of  the  probe 
sets.  The  signal  was  then  amplified  and  the  nuclei  counter- 
stained  with  1  |xg/ml  DAPI  solution  for  3  min.  Type  8  tar¬ 
get  mRNA- specific  probe  sets  with  a  550  flurophore  were 
used  as  well  as  the  Type  8  signal  amplification  systems. 
The  plate  was  imaged  and  analyzed  using  the  Spot  Detector 
Bioapplication  on  the  ArrayScan  VTI  high-content  analyzer 
(Thermofisher  Cellomics  Division,  Pittsburgh,  PA).  Using 
the  20X  objective,  the  XF93-Hoechst  and  XF93  TRITC 
filter/dichroic  combinations  were  used  to  image  DAPI  and 
the  target  mRNA,  respectively. 

SM  exposures 

SM  was  obtained  from  the  US  Army  Edgewood  Chemi¬ 
cal  Biological  Center  (Aberdeen  Proving  Ground,  MD). 
All  SM  exposures  were  performed  in  a  certified  chemical 
surety  fume  hood.  A  frozen  aliquot  of  neat  (undiluted)  SM 
in  keratinocyte  growth  medium  was  thawed  and  vortexed 
to  generate  4  mM  SM  stock  solution.  This  stock  solution 
was  placed  on  ice  and  immediately  diluted  into  medium  to 
expose  cells  to  SM.  The  exposure  dose  was  200  pM  SM 
for  all  experiments.  This  dose  was  selected  because  it  has 
been  historically  established  to  approximate  a  vesicating 
dose  in  skin  based  on  a  histopathological  comparison  of 
SM-exposed  skin  and  cultured  cells.  For  SM  exposures, 
cells  were  left  at  37  °C  with  room  C02  concentrations  for  a 
maximum  of  30  min  and  returned  to  37  °C  with  5  %  C02. 
Unexposed  controls  were  also  left  at  37  °C  with  room  C02 
concentrations  at  the  same  time.  Cell  viability  and  expres¬ 
sion  of  cytokines  were  assessed  24  h  after  exposure. 

Cell  viability  assay 

3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium 
bromide  (MTT)  was  used  to  assess  cell  viability.  A  stock 


solution  of  MTT  (Sigma- Aldrich,  St.  Louis,  MO)  was  made 
in  PBS  at  5  mg/ml,  filter  sterilized,  and  stored  at  —20  °C 
until  use.  This  stock  solution  was  diluted  to  0.5  mg/ml 
in  IX  PBS  containing  1  mg/ml  glucose  and  pre-warmed 
to  37  °C.  Twenty-four  h  after  SM  exposure,  the  cell  cul¬ 
ture  medium  from  each  well  was  removed  and  transferred 
to  96- well  U-bottom  plates,  sealed,  and  stored  at  —80  °C 
for  subsequent  analysis.  The  cells  were  then  overlaid  with 
100  pi  per  well  of  the  PBS/glucose/MTT  solution  and  incu¬ 
bated  at  37  °C  with  5  %  C02  for  2  h.  After  incubation,  the 
PBS/glucose/MTT  was  removed  and  discarded.  The  result¬ 
ant  purple  formazan,  formed  by  mitochondrial  dehydroge¬ 
nase  conversion  of  MTT,  was  dissolved  by  the  addition  of 
100  pi  per  well  isopropanol.  The  plate  was  gently  shaken 
on  a  plate  shaker  for  1  min,  and  absorbance  at  570  and 
690  nm  was  measured  on  a  SpectraMax  M5  plate  reader 
(Molecular  Devices,  Sunnyvale,  CA). 

Cytokine  analysis 

The  levels  of  IL-8  and  TNFa  in  culture  medium  samples 
from  primary  screening  and  deconvolution  were  assessed 
using  standard  capture  ELISA.  Antibodies  used  for  the 
IL-8  ELISA  were  mouse  anti-human  IL-8  monoclonal 
IgGl  clone  6217  (R&D  Minneapolis,  MN,  cat#  MAB208) 
for  capture  and  goat  anti-human  IL-8  polyclonal  (Sigma, 
St.  Louis,  MO  cat#  18026)  for  detection.  Antibodies  used 
for  the  TNFa  ELISA  were  mouse  anti-human  TNFa  mon¬ 
oclonal  IgGl  clone  28401  (R&D  cat#  MAB610)  for  cap¬ 
ture  and  goat  anti-human  TNFa  polyclonal  (R&D  cat# 
AF-210-NA)  for  detection.  Both  assays  used  the  same  rab¬ 
bit  anti-goat  IgG  (whole  molecule)  peroxidase-conjugated 
antibody  (Sigma  cat#  A4174)  for  labeling.  Capture  anti¬ 
body  was  diluted  in  PBS,  added  to  Grenier  high-binding 
white  96-well  plates,  and  incubated  overnight  at  4  °C.  The 
plates  were  then  blocked  for  2  h  at  room  temperature  with  a 
50:50  mixture  of  SuperBlock  (Thermo  Scientific)  and  5  % 
BSA  PBS  with  0.5  %  Tween  20.  The  plates  were  washed 
five  times,  samples  were  added,  and  then,  the  plates  were 
incubated  for  1  h  at  room  temperature.  Detection  and  labe¬ 
ling  antibodies  diluted  in  blocking  buffer  were  added  after 
washing  five  times  and  incubated  for  1  h  each  at  room 
temperature.  Chemiluminescent  Peroxidase  Substrate  for 
ELISA  (Sigma)  was  added  after  a  final  wash  step,  and  sig¬ 
nal  was  detected  using  a  Spectramax  M5  (auto  PMT  and 
500  ms  integration  time). 

For  target  validation,  multiplex  analysis  of  IL-6,  IL-8, 
and  TNFa  in  culture  medium  samples  was  performed  using 
a  bead-based  multiplex  cytokine  assay  (Milliplex,  Milli- 
pore,  Billerica,  MA)  and  analyzed  using  a  Bio-Plex  System 
array  reader  with  Bio-Plex  Manager  4.0  software  (Bio-Rad 
Laboratories,  Hercules,  CA).  The  data  were  expressed  as 
the  average  of  three  replicates  =b  standard  deviation  (SD) 
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over  time  ( n  =  3).  Data  were  normalized  to  the  exposed 
p38a  control  group  averaged  across  all  validation  assays 
and  expressed  as  fold  expression  relative  to  this  con¬ 
trol  group.  Data  were  analyzed  for  statistical  significance 
using  a  one-way  analysis  of  variance  (ANOVA)  followed 
by  Dunnett’s  multiple  comparison  tests  comparing  the 
exposed  negative  siRNA  pool  control  group  that  was  aver¬ 
aged  across  all  validation  assays  to  all  other  groups.  Statis¬ 
tical  significance  was  defined  as  p  <  0.05  for  all  tests. 

Results 

Primary  screening 

Prior  to  screening,  transfection  optimization  studies  were 
performed  using  siRNA  targeting  cyclophilin  b  to  establish 
transfection  conditions  necessary  to  achieve  at  least  70  % 
knockdown.  Cyclophilin  b  mRNA  levels  were  measured  in 
situ  using  the  Quantigene  View  assay  (Affymetrix,  Santa 
Clara,  CA),  and  the  results  quantified  by  high-content  anal¬ 
ysis  (HCA)  using  the  ArrayScan  VTI  high-content  analyzer 
(Thermofisher,  Cellomics  Division,  Pittsburgh,  PA).  We 
were  able  to  achieve  upwards  of  90  %  knockdown  of  cyclo¬ 
philin  b  with  as  little  as  1.2  pmol  of  siRNA  per  well  (data 
not  shown).  However,  we  elected  to  use  10  pmol  per  well  for 
the  primary  screening  and  library  deconvolution  in  order  to 
capture  potential  targets  that  might  require  higher  levels  of 
siRNA  for  effective  knockdown.  Target  knockdown  kinetics 
were  also  evaluated,  and  target  knockdown  was  maximal  at 
2  days  after  transfection  (data  not  shown).  All  transfections 
were  performed  using  plate  pipetting  robotics,  and  consist¬ 
ency  of  knockdown  across  the  plate  was  validated.  Histori¬ 
cal  microarray  data  from  SM  exposures  of  normal  human 
epidermal  keratinocytes  (NHEK)  and  the  mouse  ear  vesi¬ 
cant  model  (MEVM)  were  analyzed  using  Ingenuity  Path¬ 
way  Analysis  (IPA)  for  molecules  and  pathways  implicated 
in  the  inflammatory  response  to  SM.  This  analysis  identi¬ 
fied  2,017  molecules  which  were  selected  as  targets  for 
our  siRNA  library  for  primary  screening.  This  library  was 
screened  in  HaCaT  cells,  and  the  endpoints  assessed  were 
the  IL-8  levels  by  capture  ELISA  and  cell  viability  by  MTT 
assay.  A  positive  control,  siRNA  targeting  p38a  (MAPK13), 
was  included  on  every  plate  and  was  used  for  data  normali¬ 
zation  across  the  study.  Though  not  a  focus  of  the  current 
screening  effort,  the  effect  of  siRNA  on  cell  viability  was 
also  assessed,  and  the  results  are  shown  in  Fig.  1.  The  range 
of  the  cell  viabilities  for  all  SM-exposed  samples  was  0.95- 
52.3  %,  and  the  average  cell  viability  for  all  SM-exposed 
samples  was  8.71  %  (SD  =  6.41).  Interestingly,  six  targets 
improved  cell  viability  to  greater  than  35  %.  These  tar¬ 
gets  are  CACNA1B  (52.3  %),  ALOX5  (40.9  %),  CD2AP 
(40.4  %),  CAMK1D  (39.1  %),  CD247  (37.5  %),  and  CCR4 
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Fig.  1  siRNA  Effect  on  Cell  Viability  of  SM-Exposed  HaCaT  Cells. 
Cells  were  transfected  with  siRNA  for  48  h  and  then  exposed  to 
200  |xM  SM.  Cell  viability  was  assessed  by  MTT  assay  24  h  after 
exposure.  Six  targets  are  circled  (CACNA1B,  ALOX5,  CD2AP, 
CAMK1D,  CD247,  and  CCR4)  that  show  improved  cell  viability 
fourfold  to  sixfold  over  the  library  average  of  8.71  %.  Results  shown 
are  the  average  of  three  replicates  ( n  =  3) 
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Fig.  2  siRNA  Inhibition  of  IL-8  Production  by  SM-Exposed  HaCaT 
Cells.  Cells  were  transfected  with  siRNA  for  48  h  and  then  exposed 
to  200  |xM  SM.  IL-8  in  cell  culture  supernatant  was  assessed  by  cap¬ 
ture  ELISA24  h  after  exposure.  Shaded  area  shows  hits  within  the 
selection  metric  of  1.5-fold  the  positive  control.  Results  shown  are 
the  average  of  three  replicates  ( n  =  3) 

(37.2  %)  (collectively,  the  SD  of  these  six  targets  is  12.3). 
There  is  a  6  %  gap  in  cell  viability  between  CCR4  and  the 
next  gene,  LOC441108;  no  other  performance  gap  of  more 
than  1  %  exists  between  the  remaining  targets. 

Figure  2  shows  the  p38a-normalized  IL-8  production 
for  each  target  in  the  primary  screen.  A  hit  selection  metric 
of  1.5-fold  the  level  of  IL-8  in  SM-exposed  p38a-positive 
controls  was  chosen  as  the  threshold  for  selection  of  tar¬ 
get  for  deconvolution,  the  next  phase  of  screening.  There 
were  150  targets  that  met  this  threshold.  The  150  targets 
selected  for  deconvolution  were  mapped  using  IPA.  The 
significance  of  the  association  between  the  data  set  and 


5)  Springer 


Arch  Toxicol  (2016)  90:375-383 


379 


so  H 


■=  40  - 


20  - 


io  H 


0  J 


fr 


U) 

C7J 

o 


"D 

O 

b 

o 

*3 

d) 

c 

d> 

O 


03 
£ fl 
(5 
Oi 
tfl 

b 

fr 

cd 

E 

E 

(a 


D) 

O 

O 

-£= 

Q. 

3- 

o 


cd 

o 

Ui 

O 

Q 


di 

ui 

c 

o 

c. 

</J 

di 

□£ 

fr 

3 

CG 

E 

E 

<TJ 


0? 

E 

Q. 

O 

OJ 

> 

to 

Q 


cu 

O 


£ 

<u 

■M 

trt 

>* 

CO 

o 

O) 

_o 

& 

Id 

E 

03 

I 


03 

~o 

Jkm 

O 

□ 

<D 

m 


OJ 

> 

4-1 

L3 

03 

£ 

C 

o 

o 


[ 


£ 

to 

T3 

i— 

O 

b 


o 

(fl 


(0 

+4 

0J 

CO 


Fig.  3  Top  Ten  Canonical  Pathways  by  Ingenuity  Pathway  Analysis 


10,0 


7.5  - 


5.0  - 


2.5  - 


0,0  - 


1 1 1 1 


1 1  ■  ■ 


o 

% 

3 

pa 

E 

~to 

O 

u 


*D 

C 

C3 

a 


03 

U 

c 

w 

O 

M— 

o 

U3 

£ 

'£ 

Cd 

_C 

a 

03 

2 

i_ 

Jd 

3 

U 

a> 

Q 

m 


c 


5 

03 

to 

JO 

W 

O 

% 

O 

c 

Cm 

L3 

< 


PI 


a> 

to 

LT3 

£u 


a 

I 


nj 

0_ 


.S5 

oo 

CL 

S 

CO 


P) 

c 

c 

P) 

CO 

oo 

i 

-1 


P> 

C 


a) 

<55 

6 

03 

oc 

T3 

o 

o 

« 

o 

u 

O 

o 

_3 

a 


pj 

c 

75 

c 

pi 

<o 

c 

Id 

£ 

P 


o 

c 

7a 

c 

e> 

<55 


co 

9 

s 

o 


Fig.  4  Top  Ten  Biological  Pathways  by  Ingenuity  Pathway  Analysis 


the  pathway  was  measured  in  two  ways:  (1)  A  ratio  of  the 
number  of  genes  from  the  data  set  that  map  to  the  pathway 
divided  by  the  total  number  of  genes  that  map  to  the  canon¬ 
ical  pathway  is  displayed  (Fig.  3):  (2)  Fisher’s  exact  test 
was  used  to  calculate  a  p  value  determining  the  probability 
that  the  association  between  the  genes  in  the  data  set  and 
the  canonical  pathway  is  explained  by  chance  alone.  A  neg¬ 
ative  log  value  of  1.3  [1.3  =  -log  (0.05)]  was  set  as  being 
statistically  significant.  The  top  ten  canonical  pathways  are 
shown  in  Fig.  3.  Most  of  the  canonical  pathways  identified 
are  related  to  inflammation,  which  is  to  be  expected  given 
that  the  library  was  focused  on  inflammation  by  IPA  analy¬ 
sis  of  microarray  data.  The  top  three  biological  pathways, 
gastrointestinal  disease,  genetic  disorder,  and  inflammatory 
disease,  all  share  the  same  two  disorders,  inflammatory 
bowel  disease  (IBD)  [-log  (p  value)  =  64.1]  and  Crohn’s 
disease  [-log  (p  value)  =  59.4]  (Fig.  4). 

Deconvolution 

Because  the  primary  screening  phase  utilized  four  pooled 
siRNA  sequences  to  knockdown  each  target,  hits  from  the 
primary  screen  underwent  deconvolution  in  which  each  of  the 
siRNA  sequences  was  evaluated  separately.  As  stated  previ¬ 
ously,  150  targets  met  the  hit  selection  threshold  in  primary 
screening,  and  72  targets  have  been  deconvoluted  thus  far. 
Deconvolution  was  performed  with  HaCaT  cells,  and  the  end 
points  assessed  by  capture  ELISA  were  IL-8  and  TNFa.  Only 


two  cytokines  were  selected  for  end  points  in  deconvolution 
in  order  to  simplify  the  screening  process.  IL-8  and  TNFa 
were  selected  over  IL-6  because  they  are  more  broadly  impli¬ 
cated  in  overall  cornea  health,  whereas  IL-6  has  been  more 
frequently  implicated  in  a  specific  cornea  pathology,  neovas¬ 
cularization.  Cell  viability  was  assessed  by  MTT  assay.  All 
the  same  controls  used  for  primary  screening  were  used  for 
deconvolution.  Cytokine  data  were  normalized  to  the  p38a- 
positive  control,  and  results  are  shown  in  Fig.  5.  The  normal¬ 
ized  IL-8  expression  for  most  of  the  sequences  (323/600)  was 
less  than  1.5,  and  for  TNFa,  many  (298/600)  were  <1.  There 
were  several  considerations  for  hit  selection  in  deconvolution. 
A  primary  factor  was  that  three  or  four  out  of  four  sequences 
demonstrated  cytokine  knockdown  less  than  or  equal  to  one 
for  both  the  IL-8  and  the  combined  scores.  Other  factors 
considered  were  known  roles  for  the  target  and  likelihood 
that  the  target  was  “draggable.”  For  example,  little  is  known 
about  the  roles  of  CORIN  and  ZNF547  in  cellular  physiol¬ 
ogy,  and  these  were  selected  for  further  study.  Much  is  known 
about  the  roles  of  MAPKs  and  RELA  in  cells  and  attempts  to 
develop  drugs  for  some  of  them  have  been  largely  unsuccess¬ 
ful;  therefore,  these  were  excluded  from  study. 

siRNA  target  validation 

Twelve  targets  were  selected  from  deconvolution  for  in 
vitro  validation:  BAZB1,  CLSTN2,  CORIN,  DMBT1, 
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Fig.  5  siRNA  Inhibition  of  IL-8  and  TNFa  Production  by  SM- 
Exposed  HaCaT  Cells  in  Library  Deconvolution.  Cells  were  trans¬ 
fected  with  siRNA  for  48  h  and  then  exposed  to  200  p,M  SM.  IL-8 
and  TNLa  in  cell  culture  supernatant  were  assessed  by  capture 

GALNT10,  GFRA1,  GNAS,  HCN1,  HLA-DQA1,  IFNAR1, 
TRAF6,  and  ZNF547.  The  two  sequences  for  each  target 
with  the  best  inhibition  of  cytokine  production  from  decon¬ 
volution  were  selected  for  use  in  validation  studies.  Trans¬ 
fection  was  optimized  for  each  siRNA  sequence  studied 
using  the  Quantigene  View  assay  as  previously  described. 
End  points  were  cell  viability  by  MTT  and  multiplex  assess¬ 
ment  of  IL-6,  IL-8,  and  TNFa  by  Milliplex  bead-based 
cytokine  assay.  The  two  siRNA  sequences  for  each  target 
were  tested  separately  and  in  combination.  When  tested  in 
combination,  the  total  amount  of  siRNA  was  kept  equal  to 
that  of  separately  tested  sequences.  Experimental  groups  for 
which  both  sequences  and  the  combination  of  sequences  sig¬ 
nificantly  inhibited  cytokine  production  are  marked  with  an 
asterisk  (Fig.  6).  There  were  several  experimental  groups  for 
IL-6  and  IL-8  that  had  statistically  significant  cytokine  inhi¬ 
bition  relative  to  the  exposed  negative  pool  control  group. 
The  actual  decrease  in  cytokine  production  was  marginal, 
but  the  data  were  significant  owing  to  the  small  standard 
error.  Only  CORIN  significantly  inhibited  cytokine  produc¬ 
tion  to  less  than  half  of  the  negative  pool  control  group  for 


ELISA  assay  24  h  after  exposure.  Results  shown  are  the  average  of 
three  replicates  d=  SE  (n  =  3).  Dotted  line  shows  performance  of  the 
p38a  positive  control 


all  three  cytokines  evaluated.  For  IL-6,  GNAS  was  also  able 
to  reduce  cytokine  production  to  less  than  half  of  the  nega¬ 
tive  pool  control  group. 

Discussion 

Herein,  we  report  on  work  thus  far  on  a  high-throughput 
siRNA  screening  paradigm  to  understand  the  molecular 
mechanisms  of  vesicant  injury  and  identify  potential  tar¬ 
gets  for  the  treatment  of  this  injury.  The  siRNA  library  we 
utilized  was  focused  on  inflammatory  pathways  because 
chronic  inflammation  after  exposure  appears  to  play  a  role 
in  progressive  clinical  pathology  and  intervention  may 
improve  clinical  outcome.  The  end  points  selected  for  this 
screen,  inflammatory  cytokine  production  and  cell  viability, 
were  selected  because  targets  that  mitigate  cellular  injury 
should  reduce  SM-induced  cell  death  and  inflammation. 

Our  results  demonstrate  that  knockdown  of  GNAS 
significantly  inhibited  SM-induced  IL-6  production  and 
knockdown  of  CORIN  significantly  inhibited  SM-induced 
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Fig.  6  siRNA  Inhibition  of 
IL-6,  IL-8,  and  TNFa  Produc¬ 
tion  by  SM-Exposed  HaCaT 
Cells  in  Target  Validation.  Cells 
were  transfected  with  siRNA 
for  48  h,  and  then  exposed 
to  200  |xM  SM.  Two  siRNA 
sequences  for  each  target 
were  tested  under  sequence- 
specific  optimized  transfection 
conditions.  Each  sequence  was 
tested  separately  (labeled  as  —  1 
and  —2)  and  in  combination 
(labeled  as  —  C).  Cytokines  in 
cell  culture  supernatant  were 
assessed  24  h  after  exposure 
by  multiplex  assay.  Results 
shown  are  the  average  of  three 
replicates  =L  SE  (n  =  3).  The 
results  for  the  control  groups 
Naive,  SM,  SM  Neg  Pool,  and 
SM  p38a  are  provided  for  refer¬ 
ence.  Targets  for  which  both 
sequences  and  the  combina¬ 
tion  of  sequences  had  cytokine 
production  significantly  less 
than  the  SM-exposed  negative 
pool  control  group  are  indicated 
by  asterisks  (*p  <  0.05). 

CORIN  and  GNAS  ( underlined ) 
inhibited  cytokine  production  to 
less  than  half  the  SM-exposed 
negative  pool  control  group 
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IL-6,  IL-8,  and  TNFa  production.  The  finding  that  CORIN 
knockdown  inhibited  production  of  all  three  cytokines,  but 
that  GNAS  knockdown  inhibited  production  of  only  IL-6 
suggests  that  GNAS  may  be  downstream  of  CORIN  in 
the  regulation  of  these  cytokines  in  SM  cutaneous  injury. 
GNAS  is  a  complex  locus  which  produces  many  tran¬ 
scripts.  Gs-alpha  is  perhaps  the  most  well-characterized 
transcript,  and  it  encodes  the  alpha  subunit  of  stimula¬ 
tory  guanine  nucleotide-binding  proteins  (G  protein). 
G  proteins  play  a  wide  variety  of  physiologic  roles  and 
are  expressed  in  almost  all  cells  and  tissues.  With  regard 
to  inflammation,  Gal 2  has  been  shown  to  regulate  both 
IL-6  and  IL-8  in  oral  squamous  cell  carcinoma  (Jian  et  al. 
2013).  CORIN  is  a  serine  protease  that  converts  pro-atrial 
natriuretic  peptide  (pro-ANP)  to  biologically  active  ANP, 
which  plays  a  role  in  the  regulation  of  blood  pressure  (Wu 
et  al.  2009).  Little  is  known  regarding  other  potential  roles 
for  this  molecule.  However,  recently,  it  has  been  shown  to 
be  up-regulated  in  the  deciduas  of  the  pregnant  uterus,  and 
CORIN  gene  mutations  have  been  associated  with  preec¬ 
lampsia,  suggesting  a  role  for  CORIN  in  pregnancy  (Zhou 
and  Wu  2013).  Although  our  identification  of  a  role  for 
CORIN  in  SM-induced  cutaneous  inflammation  is  distinc¬ 
tive  from  any  previously  identified  role  for  this  molecule, 
it  is  consistent  with  its  biochemical  activity,  as  there  are  a 
number  of  serine  proteases  that  play  a  role  in  regulating 
the  immune  system  and  inflammation  (Safavi  and  Rostami 
2012). 

Although  our  siRNA  library  was  focused  on  inflamma¬ 
tory  pathways,  we  identified  six  targets  that  increased  cell 
viability  fourfold  to  sixfold  over  the  average  cell  death  for 
the  entire  screen.  Knockdown  of  CACNA1B  (the  alpha- 
1B  subunit  of  calcium  channel)  or  CAMK1D  (calcium/ 
calmodulin-dependent  protein  kinase  1 -delta)  improved 
cell  viability  following  SM  exposure.  These  findings  are 
consistent  with  the  known  role  of  calcium  signaling  in 
SM-induced  cell  death  and  differentiation  (Rosenthal  et  al. 
2000).  ALOX5  is  a  lipoxygenase  that  has  a  role  in  the 
catalysis  of  leukotrienes.  Although  prostaglandin  and  leu- 
kotriene  synthases  have  been  shown  to  have  a  role  in  SM 
inflammation  (Black  et  al.  2010;  Tanaka  et  al.  1997),  this 
is  the  first  data  to  demonstrate  a  potential  role  of  leukot¬ 
rienes  in  SM-induced  cell  death.  Interestingly,  three  of  the 
knockdown  targets  that  improved  cell  viability  have  known 
roles  in  T-cell  physiology:  (1)  CD247,  the  zeta  chain  of 
the  T-cell  antigen  receptor;  (2)  CCR4,  a  CC  chemokine 
receptor;  and  (3)  CD2AP,  a  multifunctional  adapter-type 
scaffolding  protein  involved  in  the  regulation  of  the  actin 
cytoskeleton.  This  finding  is  curious  given  that  our  studies 
were  conducted  with  HaCaT  cells,  which  are  not  known  to 
have  any  contaminating  immune  cells  or  to  express  T-cell- 
related  molecules.  It  is  important  to  note  that  these  were 
initial  findings  from  primary  screening,  and  the  role  of 


these  targets  in  SM  cell  viability  has  not  been  validated.  It 
is  also  important  to  restate  that  primary  screening  was  per¬ 
formed  using  Dharmacon  siGENOME  siRNA,  which  can 
have  greater  off- target  effects  relative  to  Dharmacon  ON- 
TARGETplus,  which  was  used  for  target  validation.  Fur¬ 
ther  studies  are  necessary  to  validate  the  roles  of  these  mol¬ 
ecules  in  SM-induced  cell  death,  and  in  the  case  of  CD247, 
CCR4,  and  CD2AP,  future  studies  should  confirm  if  in 
fact  HaCaT  cells  and  primary  keratinocytes  express  these 
proteins. 

The  top  ten  canonical  pathways  identified  by  pathway 
mapping  were  largely  related  to  immune  function,  which 
is  to  be  expected  given  the  focus  of  our  siRNA  library 
on  inflammation.  It  is  interesting  to  note  that  p53  signal¬ 
ing  and  molecular  mechanisms  of  cancer  were  two  of  the 
top  ten  pathways.  While  this  is  consistent  with  the  radi- 
omimetic  properties  of  SM,  any  potential  role  for  p53 
or  DNA  damage  in  SM-induced  inflammation  remains 
unclear  at  this  time.  The  top  three  biological  pathways  are 
all  related  by  IBD  and  Crohn’s  disease.  TNFa  is  a  major 
regulator  of  these  disorders,  and  TNFa  was  shown  to  be 
the  most  highly  interconnected  molecule  in  the  pathway 
mapping  analysis  of  our  deconvolution  library  (data  not 
shown).  This  suggests  a  similarity  between  SM  cutane¬ 
ous  injury  and  IB  D/Crohn’s  disease  at  least  in  regard  to 
an  important  role  for  TNFa  in  SM  injury.  TNFa  block¬ 
ade  is  an  effective  clinical  strategy  for  the  management  of 
IBD  and  Crohn’s  disease.  Although  TNFa  blockade  alone 
is  not  effective  in  treating  SM  cutaneous  injury,  Etaner- 
cept  (a  soluble  TNFa  receptor)  is  a  necessary  component 
of  a  combinatorial  therapy  that  has  shown  some  modest 
efficacy  in  treating  this  injury  (personal  communication, 
James  Dillman,  PhD). 

We  did  not  identify  an  effect  on  IL-ip,  which  is  com¬ 
monly  up-regulated  in  cellular  injury,  because  SM-injured 
HaCaT  cells  do  not  produce  this  cytokine  in  our  laboratory. 
Future  work  will  investigate  any  effect  on  IL-ip  produc¬ 
tion  using  primary  keratinocytes,  which  do  express  IL-ip 
in  SM  injury.  Although  we  have  many  more  targets  for 
which  the  down  selection  and  validation  process  has  yet  to 
be  completed,  we  have  identified  that  CORIN  regulates  the 
cytokines  IL-6,  IL-8,  and  TNFa  in  cutaneous  SM  injury. 
A  role  for  CORIN  in  vivo  remains  to  be  validated,  and  the 
known  role  of  CORIN  in  regulating  blood  pressure  pre¬ 
sents  important  considerations  for  its  clinical  potential  as  a 
therapeutic  target  for  SM.  However,  we  believe  this  finding 
validates  our  approach  that  high-throughput  siRNA  screen¬ 
ing  can  identify  potential  therapeutic  targets  for  chemical 
injury.  Further  studies  are  necessary  to  validate  the  roles 
of  other  targets  in  SM  inflammation,  and  additional  siRNA 
libraries,  such  as  those  focused  on  cell  death  pathways, 
should  be  screened  in  future  efforts  to  investigate  other 
mechanisms. 
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